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Abstract 

We extend a previously developed technique for computing 
spin-spin critical correlators in the 2d Ising model, to the case of mul- 
tiple correlations. This enables us to derive Kadanoff-Ceva's formula 
in a simple and elegant way. We also exploit a doubling procedure in 
order to evaluate the critical exponent of the polarization operator in 
the Baxter model. Thus we provide a rigorous proof of the relation 
between different exponents, in the path-integral framework. 
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Since Schultz, Mattis and Lieb |]J showed that Onsager's solution of 
the two-dimensional (2d) Ising model could be simply explained in terms of a 
single Majorana fermion, there has been an increasing interest in the study of 
2d statistical mechanics models by means of field-theoretical methods. In the 
same vein, Luther and Peschel @ proved that the scaling regime of the eight- 
vertex (Baxter ||) model can be described in the continuum limit in terms 
of a Thirring || Lagrangian. In this way, the 2d Ising and Baxter models 
became fruitful testing grounds for new ideas and computational methods. 
In a previous work it has been shown how to evaluate 2-point correlators 
in 2d systems [[5], through a path-integral approach to bosonisation ||. In 
particular, the critical behavior of the Ising (on-line) spin-spin correlation 
function was obtained, by using a slightly modified version of the identity 
derived by Zuber and Itzykson 0: 

F 2 2 (xi,X2) =< cr(xi)cr(x 2 ) > 2 =< exp tt / dzJ (z) > (1) 

where is the Dirac fermion current which is obtained out of the original 
Majorana fields after squaring the correlator. <> means vacuum expectation 
value (v.e.v.) in a model of free massless fermion fields. 



The purpose of this note is twofold. On the one hand we extend the 
above mentioned method to compute the 2n-point correlator. Thus, we pro- 
vide an alternative derivation of Kadanoff-Ceva's formula |j| that could be 
useful when considering certain non-trivial extensions of the Ising model such 
as the off-critical |§ and the defected |L0| cases. On the other hand we adapt 
the doubling technique f|Tl[] which led to (1), in order to calculate the cor- 
relation function of the polarization operator in the Baxter model 0. This, 
in turn, allows us to provide a path-integral confirmation of the relations 
between different critical exponents (those corresponding to energy-density, 
crossover and polarization), a result previously established by Drugowich de 
Felicio and Koberle |12| in the operator framework. 



For the sake of clarity we shall begin by briefly summarizing the main 
points of the spin-spin correlator calculation. In ref.|| the line integral in 
(1) was written as 

f 2 dzJ (z) = ( d 2 xV 

J XI J 
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where is an auxiliary vector field with components: 



A (z , zi) = 8(z )9(z 1 - x x )9{x 2 - zi) 
Ai(zo, zx) = 0. 

This simple manipulation enabled us to express the squared spin-spin corre- 
lator in terms of fermionic determinants: 

2 _ det(i0 + n£) 

F *W> X *>- det^D (2) 

where the coordinate dependence in the right hand side (r.h.s.) of (0) is, of 
course, contained in J\. 

Finally, one performs a change of path-integral fermionic variables which 
is chosen so as to decouple fermions from the background field A^. It is 
interesting to note that, in this formulation, the desired 2-point function is 



just the square root of the Fujikawa Jacobian J F |L3| associated with the 
transformation in the fermionic measure: 

F 2 (x 1: x 2 ) = J F (xi,x 2 ) 1 * (3) 
As shown in [14|, this Jacobian must be computed with a gauge-invariant 



regularization prescription in order to avoid a linear divergence (this gauge 
invariance is a consequence of a symmetry in the original lattice system ]15| ) . 
This procedure then leads to the well-known power law decay of the spin-spin 
on-line function, with exponent equal to |. 

Let us now show how to extend the above depicted technique to the 
computation of the 2n-point spin correlation function at criticality. To this 



end, we follow ref.|l6| where it was shown that, after squaring the correlator, 
each pair of consecutive spin variables can be identified with an exponential 
similar to the one appearing in (|1|) (See also ref . for a very interesting 



study on the doubling procedure and the operator content of fermion fields 
in the Ising model). We can then express the squared 2n-point correlator as 



F 2n 2 (xi, 



• %2n) 



< 



2n 

n 



a (x i) >=< 



II exp 7r / dzJ (z) > (4) 



-l,odd 
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where, as before, <> in the r.h.s. means v.e.v. to be evaluated in a model 
of massless Dirac fermions. It is apparent that each line integral in (|]) can 
be cast in the form 

x i+i r 

dzJ (z) = d zJ^A^z-.x^Xi+i) 

where we have introduced the n classical singular potentials 

A (z; Xi, x i+1 ) = 5(z )9(z 1 - Xi)6(x i+ i - z x ) 

Ai(z]Xi,x i+ i) = 0. 

In order to rewrite (£|) in a more compact way we construct a new vector 
field C M as a simple superposition of A^s: 

2n-l 

Co(z) = A {z; x h x i+1 ) (5) 

i=l, odd 

C x {z) = 0. (6) 

Thus, the 2n-point function can be expressed in terms of fermionic determi- 
nants: 

2 det(i0 + *@) 
F2n ~ deU0 ' (?) 

exactly as it happens in the n = 1 case (See @), but with A^ replaced by 
The next step is to write C M in terms of scalar functions $ c and r\ c as 



Cf, = e^dvQc + <V/c- (8) 

Now we perform a decoupling change of path-integral fermionic variables 
with chiral and gauge parameters $ c and r) c , respectively: 

* = x e -" (75 * c - i??c) (10) 
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A detailed computation of the Fujikawa Jacobian Jp associated to this change 
has been given many times in the literature (See, for instance, ref. 0); here 
we just write down the final result: 



J F = exp - / d 2 x$ c n$ c 



(11) 



We then get 



F 2n (xi, x 2n ) = Jf(xi, x 2 , x 2n ) 



(12) 



Therefore we see that, in our formulation, the squared multipoint correlator 
can be identified with a fermionic Jacobian, exactly as in the 2-point case. 
At this stage one has to solve the system of differential equations for $ c and 
rj c , obtained by replacing (g) in (^j) and (^). Finally, by inserting the result 
in ([□]) and (|12"D, one obtains 



where i > j and even (odd) refers to a constraint on i+j; i,j = 1, 2, ...2n. We 
have also set an ultraviolet cutoff, which divides the coordinate differences, 
equal to 1. This formula exactly coincides with the famous Kadanoff-Ceva's 
result g. 

Let us now study the Baxter model ||, which can be considered as two 
Ising systems interacting through their spin variables (this model is related, 
through a duality transformation, to the Ashkin- Teller model |18[]). As shown 
by Luther and Peschel |2| , the scaling limit of this model is described by the 
Thirring M interaction: 



where, as before, J M is the Dirac fermionic current and the coupling constant 
A is proportional to the four-spin coupling of the original lattice model. The 
Baxter model is known to have two natural order parameters, the magneti- 
zation and the polarization < P >=< aiSi >, where <jj and Si are the spin 
operators of each Ising system. In the continuous formulation the 2-point 
correlator for the polarization operator is given by 




(13) 



int 




(14) 



< P(x)P(y) >\ = < <? X S X (TySy > x 
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where <>a means v.e.v. with respect to the fermionic model defined by (|I4|). 
For A = the above expression becomes the squared Ising correlator. This 
suggests the following identification: 

r R 

< P(0)P(R) > x =< exp 7T / dzJo(z) > A 

Jo 

The r.h.s. of the precedent equation can be computed by employing a slightly 
modified version of the method described above. Indeed, it is easy to show 
that the introduction of an auxiliary vector field through a Hubbard- 
Stratonovich identity, allows to write 

< P(0)P(R) > A = J (15) 

with 

Z = JvA^e-I det (0 + (2 A) l ' 2 fi) (16) 

and 

Z'= [ VA„ e -f d2x ^ det (i$ + (2A) 1/2 A) , (17) 

where 



Ap = e^dv® + d^r] 

<t> B = $ + -r=$ c 
v 2A 



r/B = V+ ~^fl c - 

Let us stress that, in contrast to the previous calculation of the Ising corre- 
lator, in the present case one has to consider quantum fields $ and rj whose 
dynamics plays a crucial role in the following computation. Concerning the 
classical functions $ c and T] c , they can be determined exactly as in the Ising 
case, i.e. using formulae (||), (|6]) and (§) for n = 1. 

We shall now turn to treat the fermionic determinants appearing in flT6|) and 
( |TTD by means of decoupling changes of fermionic variables, similar to the 
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one denned by equations (§) and (|IPD, but with parameters $p and tjb in 
the form: 

V]> _ e -V2A(75*B+i??B) ^ 
\]> _ y g-v / 2A(75*s-^s) 

The corresponding Jacobian is given by 

J F = exp - f d 2 x{§ + — J=$ c )n($ + — ^=$ c ) 

TT J \/zA \/zA 



Of course, this result must be used in ([16]), whereas the same expression, but 
with $ c = is to be employed in (|17D . In so doing one readily discovers that, 
due to the fact that Jp does not depend on the field 77, this field becomes 
decoupled from <£> in both Z and Z' . As the corresponding functional integrals 
over 77 coincide, they cancelled out when performing the quotient in equation 
([15]) and one then gets 

< P(0)P(R) >a=< P(0)P(R) >o < e^J^^/A*. > 

where the first factor in the r.h.s. corresponds to the doubled Ising correlator, 
whereas the second one is a v.e.v. to be evaluated for a model of free scalars 
$ with Lagrangian density given by 



c = d + -) w.* 

I 71 



As it is well-known this computation can be done by a standard shift in the 
bosonic variable $. The final result is 

< P(0)P(R) > A = (|) 2A - (18) 

where a is an ultraviolet cutoff and Ap is the critical exponent associated to 
the polarization operator, for which we get: 



A p = (19) 



1 

41+ „ 

Recalling the results for the energy-density (e) and the crossover (CV) oper- 



ators || |12f , one obtains 

4A P = A e = (Ac,.)" 1 (20) 
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which is the relation predicted by several authors |20| and first derived 
by Drugowich de Felicio and Koberle [12] in the operator framework. 



In summary, we have extended a functional approach || , previously used 
to compute 2-point functions in 2d critical systems, to the case in which mul- 
tipoint correlators are considered. In particular, we provided an alternative 
derivation of Kadanoff and Ceva's result || for the 2n-spin on-line function. 
Our contribution can be viewed as a complement to previous works based on 
operational bosonisation, where 4-point functions were explicitly calculated 
[0] ||21|| . We feel that our formulation could be more practical when consider- 
ing, for instance non-critical correlations f9j. Indeed, in this case one expects 
to have a temperature-dependent ("massive") determinant, that can be eas- 



ily handled by following the perturbative strategy of ref . [ 22 1 . The study of 



multipoint correlators in the defected Ising model [IU| can be also envisaged 
in our scheme. 

We have also computed the 2-point function describing the critical fluctua- 
tions of the Baxter polarization operator. Thus we obtained its corresponding 
critical index. This completed the path-integral proof of the relationship be- 
tween energy-density, crossover and polarization exponents, which had been 
initiated in ref. ||. 
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